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Heavy quark hadroproduction is currently being studied both at fixed-target energies 
with proton and pion beams and at colhder energies. Additionally, charm and bottom 
quark production will be copious at the future RHIC and LHC colliders where nucleus- 
nucleus collisions are expected to produce a deconfined state of quarks and gluons, the 
quark-gluon plasma. Signatures of deconfinement include lepton pair production by partons 
in the plasma and an enhancement of heavy quark production PJ^. Thus a detailed 
knowledge of heavy quark production is required to separate its production during the 
initial nucleon-nucleon collisions from that by the quark-gluon plasma. 

Previous measurements of the cc production cross section at \^ < 63 GeV suggested 
that the lowest order (Born) cross section underpredicted the data by a factor of two to three 
I^J^, called the K factor after a similar situation in Drell-Yan production. More generally, 

where Q is the produced heavy quark. The projectile and target, A and 5, can be either 
hadrons or nuclei. The next-to-leading order (NLO) corrections to the Born cross section 
have been calculated and an analogous theoretical K factor can be defined from the 
ratio of the NLO to the Born cross sections, 

[AB ^ QQ) 

~ CTBorniAB ^ QQ) ' ^ ' 

where (Tnlo is the sum of the Born and 0{as) corrections. Particularly for the "lighter" 
heavy quarks, c and 6, the NLO cross section is strongly dependent on the choice of the 
renormalization and factorization scales which determine both -ft'exp and i^th- 

In this paper we discuss the scale dependence of c and h quark production and its influence 
on the K factors. While we explore the dependence within the limits of perturbative QCD, 
we adjust the scales and the heavy quark mass, mg, to achieve ~ iQ HI keeping in 

mind that further corrections to the cross section could also be large. We use these mass 



^^C^p*-* is the K factor for the data compared to the NLO theory calculation. 
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and scale parameters to predict the cross sections and Kth for charm and bottom production 
at RHIC and LHC energies, VS = 200-500 GeV and 5.5-14 TeV respectively. Since only 
the Born contribution is often used in event generators, it is important to check that the 
distributions calculable at the Born level are a good approximation to the NLO results. We 
show that i^th is nearly constant if the renormalization and factorization scales are assumed 
to be a function of the heavy quark transverse momentum, pt- Our calculations are done 
with a Monte Carlo program developed by Nason and collaborators [|,^|lOl- 

The double differential cross section for QQ pair production by hadrons A and B is 

^^^^ -i:!"^^ ...^-(... ..)^-(... ..)^.% ''^--'-^-;;;-;-- . (3) 

Here i and j are the interacting partons and the functions Fi are the number densities 
of gluons, light quarks and antiquarks (m < mg) evaluated at momentum fraction x and 
factorization scale fip- The short-distance cross section, dij, is calculable as a perturbation 
series in as{^iR) where the strong coupling constant is evaluated at the renormalization scale 
/iij. Both and are of the order of the heavy quark mass. 
At the Born level, the cross section can be written as 

EqEq ^^^^^^ = j -^dxi dx2 C{xi, X2) 5^(a;iPi + 0:2^2 - Pq - Pq) (4) 

where ^/s, the parton-parton center-of-mass energy, is related to V^, the hadron-hadron 
center-of-mass energy, by s = X1X2S > 'irriQ. The intrinsic transverse momenta of the 
incoming partons have been neglected. The convolution of the subprocess cross sections 
with the parton number densities is contained in C{xi, X2), 

C(x,,X2) = E[^/(^i)^#(^2) + F/(xi)Ff(x2)]^ + F/(xOF/(^^^ , (5) 

where production is by qq annihilation, qq QQ, and gluon fusion, gg — > QQ. The scale 
dependence has been suppressed since there is no distinction between fip and fiR at this 
order. The scale 2mQ is commonly used for the total cross section, motivated by the s- 
channel processes. However, the t and u channel gluon fusion graphs involve heavy quark 



exchange between the gluons, suggesting that tuq is a better scale choice. Thus some 
ambiguity in scale already exists at leading order. 

Four-momentum conservation leads to the rather simple expression 

-^^ = x,x.C(x„x.), (6) 
dp^dygdyQ 

where the momentum fractions, Xi and X2, are 

xi,2 = ^(e±^'5+e±^) , (7) 

and rriT = + pj^. The target fraction, X2, decreases with rapidity while the projectile 

fraction, xi, increases. Both increase with pt- The quark px distribution is determined by 
ddij/dt oc as discussed in ||TT|. 



At NLO, 0{al), in addition to real and virtual corrections to the Born diagrams, quark- 
gluon scattering, q{q)g QQ(l{q), is also included. At the Born level, this has been in- 
terpreted as scattering a heavy quark in the hadron sea with a light quark or gluon and 
referred to as flavor excitation| [|ri|]. The relative importance of the excitation diagrams to 
hadroproduction is shown in Ref. . 



The total short- distance cross section, dij, is 



a.,{s,mQ,^F,^n) = {/"(p) + 47ra,M [/^.(p) + ^.(p) ln(p|/my] + 0{al)} , (8) 

where p = Amg/s and the functions flj are coefficients of the perturbative expansion. The 
Born contribution is given by ffj and vanishes when ij = gq, gq. The Born coefficients 
are always positive and ffj — > as p ^ 0. The NLO coefficients can be either positive or 
negative and tend to finite values as p — > in the gg and qq channels. The mass factorization 
terms, flj, only contribute when p/j 7^ ttiq. The gg and qg contributions are dominant in 
the high energy limit due to t-channel gluon exchange. The NLO gg corrections are very 



^In deep-inelastic scattering when S> mg the heavy quark can be treated as massless and 
absorbed into the parton densities but when ~ mg, flavor creation diagrams are dominant |12|. 
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large since fgg/fgg ^ oo as p — ^ 0, indicating that a small x resummation is needed. 
For completeness, in Fig. 1, the coefficients are reproduced from Ref. as a function of 
Vv^ = V^/2?TiQ. Although the coefficients are shown for < 200 only, flj and fl^ have 

reached a plateau. The maximum value of l/y/p for a given energy is 1 / y/pma.x = ^fS jlrrtQ^ 
so that 3 < 1 / ^Pmax < 12 at typical fixed-target energies while 1 / ^Pmax > 20 for 66 
production and > 77 for cc production at RHIC and LHC. 

The physical cross section should be independent of the scale: the dependence in eq. 
(8) introduces an unphysical parameter. If the perturbative expansion is convergent, i.e. 
if further higher-order corrections are small, at some scale the contribution to 

the cross section should be smaller than the contribution^]. If the scale dependence 

is strong, the perturbative expansion is untrustworthy [|ll|]. Since i^th — 1 > 1, further 
higher-order corrections are needed, particularly for charm and bottom quarks which are 
rather "light" when is large. Although the scales are, in principle, independent, we take 
fJ'F = fJ'R = fJ' unless otherwise noted because this assumption is inherent in global analyses 
of parton densities. Sometimes we use the notation /i = \/Q^. 

The two loop value of the strong coupling constant, as, depends on the number of active 
flavors, /, and the appropriate value of Aqcd for the number of flavors, Af, 

6'^lnln(/iVA2 



hf ln(/i2/A2 



(9) 



hf ln(/iVA 

where hf = (33 - 2/)/127r and h'f = (153 - 19/)/(27r(33 - 2/)). For charm and bottom 
production, / = 3 and 4. Parton densities are calculated assuming the heavy quark con- 
tributions turn on at /i = mg and as{mQ, f) = asimq, / + 1). Above this threshold, some 



groups, e.g. treat the heavy quark as massless in the evolution of the parton densities, 
introducing some overcounting into production calculations. Also note that the heavy quark 



^The order of the expansion is represented by n where n>2 for QQ production. A calculation to 
order 0{d^) introduces corrections at the order 0(a"^^). Thus the fi dependence should decrease 
when additional higher-order corrections are included if the perturbation theory converges. 
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mass threshold chosen in the global analysis of parton densities may differ from the mass 
used in calculations. 

We have used two sets of recent parton distribution functions^, GRV HO ||16| and MRS 
D— ' |T^. GRV HO has a low initial scale, Qqgrv ~ 0-3 GeV^, with valence-like parton 



distributions, therefore evolving very quickly with Q^. MRS D-' has a more conventional 
initial scale, Qqmrs ~ ^ GeV^, and sea quark and gluon densities that grow as ~ 
when X ^ 0. The difference in Ql produces the contrast in QQ production between the two 



sets. Both are compatible with the recent deep-inelastic scattering data from HERA |[17 



although both are on the high side of the data. We also use the SMRS P2 |18| and GRV HO 



pion |jT9[ densities. The GRV HO pion distributions are obtained from their proton set and 
have a similar behavior at low x. The SMRS P2 distributions are based on an older set of 
nucleon distributions with a different value of A4 than MRS D— '. In this case, we evaluate 
Os at the MRS D— ' value of A4. Also, note that at Qqmrs' pion sea quark and gluon 
distributions are assumed to be constant as x — > 0, incompatible with the MRS D— ' low x 
behavior. 

Since we will discuss both p{p)p and 7r~p production, it is instructive to show the quark 
and gluon momentum distributions for protons and pions. In Fig. 2(a) and (c) we show 
xf\x) = x{uv{x) + dv{x) + 2(u{x) + d{x) + s{x))) = x{uv{x) + dv{x) + S{x)) and xg{x) 
at = 5, 25, and 100 GeV^. The lower curve for each set corresponds to the lowest Q^, 
typical for charm production. The middle curve, at 25 GeV^, represents charm production 
at pt — 3.5mc and b production at /x ~ m^. The highest is equivalent to c and b 
production at pt > ttiq. In 2(b) and (d), we show the parton distributions as a function of 

for X = 0.0025, 0.007, and 0.15. Here the lowest x values correspond to b production at 
CDF and UAl. We choose x = 0.15 as typical of central charm production at fixed target 



^All available parton distribution functions are contained in the package PDFLIB |15|, available 
in the CERN library routines. 
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energies. The curves are MRS D-' (solid), GRV HO (dashed), SMRS P2 (dot-dashed), 
and GRV HO pion (dotted). At fixed Q^, parton distributions that initially increase as 
X are depleted at higher x compared to constant initial distributions such as SMRS 
P2. The pion quark distributions are larger than the proton distributions as a; — > 1. When 
the partonic cross sections are convoluted with the parton densities, gg and qg processes 
dominate because the fiux, oc F^^{xi, fiF)F^ {^2, fip)^, is large. In fact, although the gluon 
density is larger at small x, as x increases, it drops faster than the quark density so that, at 
some X1X2, FqFg > F^F^ . Eventually, as X1X2 1, qq annihilation dominates since the 
valence quark density is most important at large x. 

In Fig. 3(a), we show the scale dependence of the NLO calculations with the pp and 
pA data on a]^{S) [^|5|,p0| assuming a linear nuclear dependence ||21|]. The scale is varied 



between mQ/2 and 2mQ with 1.2 < niq < 1.8 GeV to show the range of theoretical uncer- 
tainty 1,13. We further assumed that a{Ds) ^ 0.2 a{D^ + D+) and a{Ac) ~ 0.3a{D° + D+) 



so that alf' ~ l.bcxlDD), as in Ref. Since rric < <5o,mrS) we take /U_f = for the 

calculations with the MRS parton densities. The three solid curves are calculated with MRS 
D— ' densities and = 1.2 GeV, /x/j = mc/2 (upper); = 1.2 GeV, /i/j = 2mc (middle); 
and rric =1.8 GeV, /i/j = 2mc (lower). The difference between the upper and lower curves 
gives the maximum variation of o"*^* for these densities, a factor of 90 for = 20 GeV 
and 20 for = 14 TeV. If the MRS DO' set is used at the same mass and scale, the cross 
section is slightly larger at fixed-target energies because the DO' gluon density is larger for 
X > 0.06 but since xg{x 0, mrs) ^ constant, a*^t(D-')/(T*§*(DO') ^ 20 at 14 TeV. 
The corresponding K factors are given by the solid line in Fig. 3(b)-i^th has very little scale 
dependence with these parton densities. 

Because the GRV HO distributions have a much lower Qq, we take fiji = fip. The dot- 
dashed and dotted curves show fi = mc/2 and fi = 2mc- The upper curves have = 1.2 



^In this notation, e.g. Fg{x) = g{x). 
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GeV, the lower, nic = 1.8 GeV. Even though the cross section is larger at low \/S for 
/i = mc/2, as the energy increases, the small x gluon density at the higher scale becomes 
the dominant feature, causing the crossover shown at high energy. In fact, because the 
gluon density starts out valence-like at Qo.grV) it is almost constant for x values probed at 
75 > 100 GeV for /i = 0.6 GeV, causing the sudden flattening of the upper dot-dashed 
curve. The variation in i^th is large when yU = mc/2. When = 2mc, -ft'th is approximately 
the same as the MRS distributions. However, the maximum variation of (T*§* is smaller for 
the GRV HO densities, a factor of 60 for = 20 GeV and 6 for = 14 TeV (excluding 
^^ = 0.6 GeV). 

Previously 0, the NLO calculations were compared to the data to fix rric and fi at 
K^J^^ ~ 1 to provide an estimate that could be extrapolated to nuclear collider energies. 
Reasonable agreement was found for = 1.2 GeV, /i = 2mc for MRS D— ' (central solid 
curve) and = 1.3 GeV, fi = m^. for GRV HO (dashed curve)^. Note however that both 
curves tend to underestimate a"^^ with K^-^ ~ 1.1 — 2. In the range of the parameter space 
defined by mg, /i/j and iip^ can be reduced to unity. However, it is questionable if the 

mass and scale values needed for ~ 1 are consistent with a perturbative treatment 

and with the defined limits of the parton density distributions^. It is also not clear that the 
NNLO corrections to heavy quark production would not be at least as large as the NLO 
corrections, particularly when mq <^ v^, even though for high-mass Drell-Yan production 

comparison to the cc data with vr" beams flpl,^ using the same parameters gives agreement 
with the data at a similar level. 
^Recently, rric = 1.5 GeV was found to be compatible with this data with some essential caveats: 
and were varied independently and out-of-date parton distributions fit with several values of 
^QCD were used . Decreasing /i/j with respect to and increasing Aqcd result in significantly 
larger cross sections for a given rrif.. Additionally, different parton densities were used in the 
calculations of c^^^y and high energy h production. 
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at NNLO ctnnlo/cnlo ~ 1.1-1.3 p2|, due to cancellations among the different channels. 

However, the variation in cr*|* is much smaller than For /i = m^, a {nib = 

4.5GeV)/(T(m6 = 5 GeV) ^ 25 at = 40 GeV and 1.4 for 14 TeV while for mb = 4.75 
GeV, a{fib = 0.5mb)/a{nb = 2mb) ^ 3 at = 40 GeV and 1.2 at 14 TeV. The K factor is 
also smaller; Kth ~ 1.2 — 2.5. 

In Fig. 4 we show the scale variation of a*^* for ir^p production at 340 GeV (a) and pp 
production at 800 GeV (b) using the GRV HO pion and nucleon distributions. As fi/nic 
increases at fixed energy, the cross section varies less rapidly but Kth increases for /i > rric. 
There is no value of fi where the NLO corrections to the Born cross section are minimal, i.e. 



no optimal scale p3|] . 

Although there is no physical reason for assuming that and /i/j are different because 
they are not separated in analyses of the parton densities, we also show the change of 
the cc cross section in tt~p production with p.„ = 500 GeV induced by varying the scales 
independently: /i^? = ttIc in 4(c) and fiR = rric in 4(d). In addition to the GRV HO results, 
we include the calculations with the MRS D— ' parton distributions and = 2mc in (c) 
and fiR = 2mc in (d) at NLO (dot-dashed) and the Born level (dotted). The running of 
is the source of the strong dependence observed in (c). When fi^ < nic, the coefficients f]j, 
usually negative, produce an overall enhancement of dij. However for fin > rric, the variation 
of Bij with fi decreases since the NLO corrections partially cancel each other. There is no 
strong effect on Kth. In (d), the increase in (T*§* for yUi? < rUc with the GRV HO distributions 
occurs because at yU^ ~ <5o,grv and low x, the sea quark and gluon distributions are valence- 
like [0. The MRS D— ' variation is negligible since fip < Qo,mrs for I^f/i^c < 2.7, thus 

IJ'F = Qo,MRS- 

Until now, we have discussed a^°^ where /i oc itlq. However, for single inclusive or double 
differential distributions, it may be more appropriate to choose a scale proportional to the 
transverse momentum of the heavy quark, pp, or its transverse mass, mp. The transverse 
momentum is a natural scale in jet and prompt photon production with massless final-state 
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particles. Additionally, since pr, ttit oc (t — mQ){u — mQ)/s — (T — mQ){U — mQ)/S, this 
choice is invariant under scale transformations of the initial state momenta and results in 
a simple relation between the partonic and hadronic variables, unlike the choice /i = s. A 
constant scale would be appropriate if mq ~ pr, but generally rric ^ pr and nib ^ Pt 
at collider energies. Therefore the pt dependent scale absorbs (resums) large logarithmic 
terms such as \n{pT/mQ) appearing when pt ^ ttlq and producing collinear divergences 
which are unregulated if /i = rriQ. We examine the pt distributions of b quark production at 
collider energies for two specific constant scale, /i oc m?,, and running scale, /i oc m^, 

to see if the data favors a particular scale. We also show the pt distributions and Kthipr) 
for charmed quarks at fixed target energies. 

We first compare the constant and running scales for single b production from the pp 
colliders at ^/S = 630 GeV and 1.8 TeV. The measurements are integrated over px above 
each pT,min- We use = 4.75 GeV. The NLO calculations with the GRV HO distributions 
are shown in Fig. 5(a) for UAl ^ and Fig. 5(b) for CDF ^ and DO |26|. The solid 



curves show ^ = rriT, the dashed curves, fi = rrif,. Neither are clearly favored by the data. 
The Tevatron data is also compared to calculations with /i = itit/^ and rrih/A. However 
H = is in clear contradiction with the data. 

In Fig. 6(a) we show the charmed quark pt distributions for Tc^p production with piab = 
500 GeV using the GRV HO distributions. The solid and dashed curves are the NLO and 
Born results for /i = rriT while the dot-dashed and dotted curves are the calculations with 
/i = rric. The NLO result drops below the Born for px > ^ GeV. High statistics charm 
measurements should be able to detect a decrease in the slope of the px distribution for 

> 4 GeV if nature favors the constant scale. No such change has been observed at 
present although data in this region is rather scarce. The resulting ratios i^'th(PT) are shown 
in Fig. 6(b) for the running scale (solid curve) and the fixed scale (dashed curve). When the 
scale increases with pt-, -ft'th is relatively constant or decreases slightly. However, the fixed 
scale produces a strong decrease of K^yi with pt- 

The calculations in Fig. 6(a) and (b) are particularly illustrative of both how the pro- 
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duction processes compete with each other and the relative importance of the NLO correc- 
tions for each process. The crossover of the Born and NLO pt distributions with the fixed 
scale, occuring at pt ~ 4 GeV, can be understood by an examination of the coefficients 
of the perturbative expansion. To schematically illustrate the pr dependence, we can take 
~ V^I'^mT in Fig. 1. In the region < < 5 GeV, 1/ ^/p < 12. For 1/^ < 4, 
f^g > fgg whUc /^^ domluatcs for l/^/p > 4. Both coefficients are always positive although 
f^gg is negative for 1/y/p > 2.2. The coefficient fgg is positive when l/-s/p > 2.2 although 
smaller than fgg-, f^g is negative and fgg > \flg\ for most values of 1/a/p- The NLO cor- 
rection to qq annihilation, /^-, is negative for 1.4 < l/s/p < 4.5. When n ^ mg, /^^ can 
partially compensate for this negative contribution. As pr increases, i/s/p enters the region 
of negative fgg at ~ 4 GeV while is large and positive. At the same time the Born 
contribution to gg fusion becomes larger than the NLO correction. These competing terms 
produce the crossover observed in Fig. 6(a) for the constant scale. The same effect causes 
the turnover in the dashed curve of Fig. 5(a) for pT,min ~ 50 GeV. Thus the constant scale, 
although introducing divergences, shows the structure of the NLO corrections. 

Assuming that // oc uit softens these effects, partially because of the f]j contributions 
but mainly due to the evolution of the parton densities. As pt increases, Xi and X2 also 
increase so that while the evolution increases the parton density at fixed x, the higher 
X values actually deplete the parton flux. This relative balance produces a nearly constant 
Kth. Results with the MRS distributions are shown in Fig. 6(c) and (d). Since fi — 2mc for 
the fixed scale, a larger pt is needed to enter the crossover region in 1/y^. 

We now turn to the scale dependence of heavy quark distributions at nuclear collider en- 
ergies. While the NLO calculations are needed for the pr dependence of QQ pair production 
and decay, trivial at the Born level, it would be convenient if the other relevant distributions 
could be modeled by the Born distributions to within a constant K factor. The QQ pair 
distributions are essential to determine the heavy quark contribution to lepton pair produc- 
tion. In Figs. 7-10 we show Kth for = 200 GeV (RHIC) and 5.5 TeV (LHC) for c and b 
production assuming /i — uit for single inclusive quark xp, rapidity, and p^ distributions in 
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(a), (c), and (e) respectively. The QQ pair xp, y, and invariant mass distributions are shown 



in (b), (d), and (f). For the pair distributions, we assume fi oc ^ttiq + {Pt,q '^p'tq)/'^- "^^^ 
results with the GRV HO parton densities are given by the circles while the diamonds are the 
MRS D— ' results. (For the NLO distributions themselves, see Ref. 0.) With the running 
scale, Kth is indeed nearly constant for these quantities, even in a regime where itlq <^ ^/S. 
Some variations occur near the boundaries of phase space, seen in the c and b quark rapidity 
distributions at y/S = 200 GeV. The K factors are independent of parton density except 
for charm production at = 5.5 TeV where Kth is larger for the GRV HO distributions. 
For the charmed quark, KthiPr) increases 50% at = 200 GeV and a factor of two at 
= 5.5 TeV. These increases perhaps indicate the appearance of large logarithms that 
need to be resummed in the pt distributions. Although some variation appears, event gen- 
erators for heavy quark production can, with relative confidence, scale all non-trivial Born 
result s0. 

Thus for heavy quark and quark pair distributions calculable at the Born level, Kth is 
nearly constant provided that fi oc my. The actual value of Kth can be determined by a 
comparison of the NLO and Born total cross sections. We have set the scale by making KfJ^'~' 
close to unity to better estimate production at nuclear colliders. Other, more sophisticated 
methods can be used to optimize the scale or relate it to other scales in perturbative QCD 



Further theoretical advances can perhaps also reduce the uncertainties. A resummation 



of the soft gluon contribution to top production near threshold has been performed |2^ . This 



near-threshold resummation may also be applied to charm and bottom ||30| production at 
fixed-target energies, suggesting a more appropriate scale. 



constant scale is assumed in [27|, leading to the conclusion that Kthipr) is a strongly increasing 
function of px- Our calculations with the constant scale confirm these results in the pT region 



explored in |27|. However, we note that for pT > 6-8 GeV at RHIC Kth decreases again while at 



the LHC Kth reaches a plateau when pT ~ 10 GeV. 
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An understanding of heavy quark production is essential to correctly interpret data from 
RHIC and the LHC. At these nuclear colliders, QQ production is in the high-energy limit 
where large logarithmic corrections to Bij{p) are important. A first attempt has been made 
to resum the small x corrections to the heavy quark production cross section in the p — > 



regime |^ which could lead to an improvement in the RHIC and LHC estimates. 
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Figure Captions 



Figure 1. The coefficients of tlie perturbative expansion, /(p), given in Ref. and shown as 
a function of The solid line is the Born contribution, /°-, the dot-dashed and dashed, 

the NLO corrections flj and f]j for (a) qq annihilation, (b) gluon fusion, and (c) quark-gluon 
scattering. 

Figure 2. The parton distributions xf{x) = x{uv{x) + dv{x) + S{x)) and xg{x) as a function 
of X and Q^. In (a) and (c), = 5, 25, and 100 GeV^ for the lower, middle, and upper 
curves respectively. In (b) and (d), x = 0.0025, 0.007, and 0.15 in the upper, middle, and 
lower curves respectively. The curves are MRS D— ' (sohd), GRV HO (dashed), SMRS P2 
(dot dashed), and GRV HO pion (dotted). 

Figure 3. The variation in cr*§*(S') and i^th('S') with parton density, rric and fi. In (a) the 
three solid curves are calculated with MRS D— ' densities and = 1.2 GeV, /i/j = mc/2 
(upper); nic = 1.2 GeV, fiR = 2mc (middle); and rric = 1.8 GeV, fiR = 2mc (lower). The 
other calculations are with the GRV HO densities. The dot-dashed and dotted curves show 
H = mc/2 and p = 2mc- The upper set has = 1.2 GeV, the lower, rric = 1-8 GeV. 
The dashed curve is = 1.3 GeV, fi = rric. -^th is shown in (b). The solid curve is MRS 
D— rric = 1-2 GeV, /i = 2mc, the dotted curve GRV HO, rric = 1-2 GeV, fi = 2mc. The 
dot-dashed curves are GRV HO, fi = O.drric, and rric = 1-2 GeV (upper), 1.8 GeV (lower). 

Figure 4. The variation of o"*§* with fi = fip = fiR for 7r~p production at 340 GeV (a) 
and pp production at 800 GeV (b) using GRV HO. The change of the cc cross section in 
TT^p production with = 500 GeV from varying the scales independently is shown for 
fip = nrriQ in (c) and for /i/j = nrnq in (d) where n = 1 for GRV HO and 2 for MRS D— '. 
The curves are GRV HO NLO (solid). Born (dashed), MRS D-' NLO (dot-dashed). Born 
(dotted). 

Figure 5. Single b production at NLO is compared with data from UA 1 (a) and CDF (b). 
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The solid curves show /i — rriT and the dashed curves, ^ — mi,. We also present // = rriT/^ 
and 777.6/4 in the upper solid and dashed curves of (b). 

Figure 6. The charmed quark pt distribution for vr^p production with piab = 500 GeV is 
shown in (a) for the GRV HO parton distributions. The sohd and dashed curves are the 
NLO and Born results for = rriT while the dot dashed and dotted curves are for /i — rric- 
In (b) we show the K^^ factors for the running scale (solid curve) and the fixed scale (dashed 
curve). The results with MRS D— ' distributions are given in (c) and (d). 

Figure 7. The ratios K^^ for charm production at \/S — 200 GeV assuming /j, = mx for 
single inclusive quark xp, rapidity, and distributions in (a), (c), and (e) respectively. The 
QQ pair xp, y, and invariant mass distributions are shown in (b), (d), and (f). The GRV 
HO results are given by the circles, MRS D— ', the diamonds. 

Figure 8. The same as Fig. 5 for charm production at y/S = 5.5 TeV. 

Figure 9. The same as Fig. 5 for bottom production at VS — 200 GeV. 

Figure 10. The same as Fig. 5 for bottom production at y/S — 5.5 TeV. 
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